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ON  THE  KINEMATICS  AND  CONTROL  OF  WHEELED  MOBILE  ROBOTS 
J.  C.  Alexander  and  J.  H.  Maddocks 


A  wheeled  mobile  robot  is  here  modelled  as  a  planar  rigid  body  that  rides 
on  an  arbitrary  number  of  wheels.  The  connections  between  the  rigid  body 
motion  of  the  robot,  and  the  steering  and  driving  controls  of  wheels  are 
developed.  In  particular,  conditions  are  obtained  that  guarantee  that  rolling 
without  skidding  or  sliding  can  occur.  Explicit  differential  equations  are 
derived  to  describe  the  rigid  body  motions  that  arise  from  rolling  trajectories. 
The  simplest  wheel  configuration  that  permits  control  of  arbitrary  rigid  body 
motions  is  determined.  The  question  of  slippage  due  to  misalignment  of  the 
wheels  is  investigated  based  on  a  physical  model  of  friction.  Examples  are 
presented  to  illustrate  the  models. 


1.  Introduction 

In  this  paper,  we  analyze  the  behavior  of  a  wheeled  mobile  robot,  or  WMR.  Such  robots  ride 
on  a  system  of  wheels  and  axles,  some  of  which  may  be  steerable  or  driven.  There  are  many 
wheel  and  axle  configurations  that  have  been  used  for  WMRs  (Balmer,  Jr.  [1982];  Carlisle  [1983]; 
Everett  [1979];  Giralt,  Sobek  &  Chatila  [1979];  Helmers  [1983a];  Helmers  [1983b];  Helmers  [1985]; 
Holland  [1985];  Hollis  [1977];  Ichikawa,  Ozaki  &  Sadakane  [1983];  Iijima,  Kanayama  &c  Yuma 
[1981a];  Iijima,  Kanayama  &  Yuma  [1981b];  Johnson  [1984];  Lewis  &  Bejczy  [1973];  Marrs  [1985]; 
Moravec  [1980];  Moravec  [1983];  Moravec  [1986];  Nilsson  [1984];  Podnar,  Dowling  &  Blackwell 
[1984];  Rogers  [1984];  Smith  &  Coles  [1973];  Wallace  et  al.  [1985];  Whitaker  [1962];  Wilson  [1985]). 
The  ultimate  objective  of  our  investigations  is  the  complete  description  of  the  kinematics  and 
the  control  theory  of  such  robots  for  low-speed  maneuvering. 

Much  of  the  research  cited  above  is  described  in  a  recent  study  of  WMRs  (Muir  &  Neuman 
[1986]).  This  work  develops  a  formalism  that  is  used  first  to  model  the  kinematics  of  each  wheel, 
and  second  to  amalgamate  the  information  about  individual  wheels  to  describe  the  kinematics 
of  the  WMR  regarded  as  a  whole.  A  condition  is  developed  that  determines  whether,  given  the 
configuration  of  the  wheels  in  the  WMR,  pure  rolling  is  possible.  If  not,  a  least-squares  fit  to 
rolling  is  obtained. 

Some  of  the  same  questions  are  considered  here.  However,  our  approach  is  somewhat  dif¬ 
ferent,  and,  we  believe,  complementary.  Attention  is  restricted  to  the  problem  of  maneuvering 
a  WMR  on  a  horizontal  plane.  Precise  and  explicit  connections  between  the  steering  and  drive 
rates  of  the  various  wheels,  and  the  position  and  orientation  of  the  robot  are  obtained.  The  in¬ 
verse,  or  control,  problem  of  determining  the  steering  and  drive  rates  that  produce  a  prescribed 
robot  trajectory  is  also  resolved.  We  determine  the  simplest  configuration  of  steerable  and  driven 
wheels  that  allows  the  robot  to  be  maneuvered  in  arbitrary  planar  motions.  Slippage  due  to  wheel 
configuration  is  also  considered.  Because  our  description  of  slippage  is  based  on  mechanics,  it  is 
necessary  to  consider  friction,  for  which  the  simplest  model  is  Coulomb’s  Law.  This  law  leads  to 
the  minimization  of  a  nonsmooth  friction,  or  dissipation,  functional  which  is  the  sum  of  absolute 
values,  rather  than  a  least-squares  functional  as  considered  by  Muir  &;  Neuman  [1986],  The 
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friction  functional  is  not  well-understood,  and  in  this  paper  a  complete  analysis  is  presented  for 
only  the  simplest  cases. 

The  qualification  of  low  speed  arises  in  this  work  because  our  model  of  rolling  does  not 
explicitly  balance  forces  and  accelerations.  It  will  be  apparent  that  this  rolling  model  retains  its 
validity  until  the  inertial  forces  arising  from  accelerations  are  so  large  as  to  saturate  the  available 
frictional  forces  between  the  wheels  and  the  surface.  Thus  the  low-speed  model  may  in  fact  be 
valid  for  relatively  fast  motions,  provided  either  that  the  turns  are  not  too  tight,  or  the  friction 
between  the  wheels  and  the  surface  is  sufficiently  large. 

In  §2,  notation  and  definitions  are  introduced,  and  some  basic  general  results  are  presented. 
In  the  first  instance,  all  wheels  are  assumed  to  roll;  i.  e.,  there  is  no  slipping  or  sliding.  This 
requirement  places  compatibility  conditions  on  the  motions  of  the  various  wheels.  These  con¬ 
ditions  and  their  consequences  are  discussed  in  §3,  which  is  a  reformulation  and  extension  of 
results  obtained  in  Alexander  &  Maddocks  [l988j.  In  this  prior  work,  attention  was  focussed 
on  intrinsic  properties  of  rigid-body  trajectories,  such  as  curvatures  and  centers  of  rotation.  In 
the  present  development,  more  emphasis  is  given  to  formulations  that  allow  efficient  numerical 
treatment.  The  conclusion  of  §3  is  that  —  disregarding  issues  of  mechanical  stability  —  two 
wheels  that  are  both  driven  and  steered  are  necessary  and  sufficient  for  the  control  of  arbitrary 
planar  motions.  On  the  other  hand,  the  four  controls  (two  steering,  two  driving)  must  satisfy 
one  (transcendental)  compatibility  condition.  The  explicit  connections  between  the  controls  and 
the  rigid-body  motion  are  described  in  §4,  where  a  resolution  of  the  inverse  problem  is  also  given. 
A  practical  design  is  suggested  in  §5,  and,  in  §6,  the  kinematics  of  three  existing  WMRs  are 
compared  with  our  proposed  design. 

The  kinematic  development  is  then  extended  to  consider  failure  of  the  rolling  model  due 
to  slippage  of  the  wheels.  There  are  two  distinct  circumstances  in  which  slippage  will  occur. 
The  first  has  already  been  mentioned  —  the  rolling  model  can  fail  to  be  a  good  approximation 
because  of  large  inertial  forces  that  saturate  the  available  friction.  This  mode  of  failure  of  the 
rolling  model  is  associated  with  high  speed  maneuvering.  The  second  mode  of  failure  is  that 
the  steering  and  driving  controls  of  the  wheels  are  not  compatible  with  pure  rolling.  This  lack 
of  compatibility  can  arise  at  any  speed.  The  first  circumstance,  which  we  call  skidding,  is  not 
considered  in  this  paper.  However,  quasi-static  evolution  equations  for  low  speed  maneuvering 
involving  the  second  mode  of  failure,  which  we  call  slippage  (sometimes  called  scrubbing),  are 
developed  in  §7  and  illustrated  in  §8. 
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2.  Notation  and  definitions,  elementary  results 

A  WMR  is  modelled  as  a  planar  rigid  robot  body  that  moves  over  a  horizontal  reference  plane 
on  wheels  that  are  connected  to  the  body  by  axles.  We  denote  the  unit  vertical  normal  by  k. 
For  our  purposes,  the  only  role  of  the  body  of  the  WMR  is  to  carry  a  moving  coordinate  system, 
the  body  coordinates  x  =  (xi,  x2),  of  a  point.  The  underlying  two-dimensional  space  coordinates 

are  denoted  X.  The  ith  axle ,  i  —  1, _ ,  m,  is  attached  to  the  body  at  the  axle  or  constraint 

point  x,-.  The  axle  is  supported  by  a  single  simple  wheel,  which  is  idealized  as  a  disc  without 
thickness  of  radius  R,  that  lies  in  a  vertical  plane  through  the  axle  point.  A  wheel  can  rotate 
in  its  vertical  plane  about  its  center  point  (which  is  attached  to  the  axle).  If  it  is  driven,  it  will 
rotate  at  a  prescribed  angular  speed.  Otherwise  it  is  passive,  and  its  rotation  rate  is  determined 
by  the  kinematics  of  the  WMR.  It  may  also  be  possible  for  the  axle  to  rotate  in  the  WMR 
about  the  vertical  through  the  axle  point.  If  the  wheel  (or  axle)  is  steered,  this  rotation  rate  is 
prescribed.  Otherwise  the  wheel  is  unsteered.  A  fixed  wheel  is  one  for  which  the  axle  cannot 
rotate.  Technically,  a  fixed  wheel  is  a  special  type  of  steered  wheel,  but  it  is  useful  to  maintain 
a  distinction. 

There  are  a  number  of  other  types  of  wheels  which  are  adopted  in  robot  design,  of  which 
a  simple  wheel  (called  a  conventional  wheel  in  §3  of  Muir  &  Neuman  [1986])  is  but  the  most 
straightforward.  Our  kinematic  analysis  can  be  expanded  to  include  many  of  these  more  com¬ 
plicated  wheels.  Some  remarks  along  this  line  are  included  in  §§5  and  6,  where  practical  wheel 
configurations  for  WMRs  are  considered. 

Regarded  as  a  rigid  body  in  three-dimensional  space,  a  simple  wheel  has  a  vector-valued 
angular  velocity  Since  the  wheel  remains  in  a  vertical  plane  at  all  times,  the  horizontal 
component  of  u,,  denoted  a<  =  o>,  —  (k  •  «<)k,  is  perpendicular  to  the  plane  of  the  wheel.  The 
vector  a,-  is  called  the  axle  vector.  The  magnitude  pi  of  the  axle  vector  is  the  rotation  rate  of 
the  wheel.  The  axle  vector  a,-  makes  a  steering  angle  0,-  from  a  reference  direction  fixed  in  the 
body  of  the  WMR,  and  an  angle  0;  from  a  reference  direction  fixed  in  space.  The  angles  0,-  and 
©,  are  well-defined  whenever  pt-  ^  0.  The  conventions  adopted  here  assume  0,-  to  be  interpreted 
modulo  2x  and  pi  to  be  nonnegative.  In  practical  implementations,  it  may  be  more  convenient  to 
interpret  d{  modulo  ir  and  to  allow  pi  to  be  negative  (i.e.,  the  wheel  may  rotate  backwards).  We 
remark  that  the  definitions  of  0<  and  0<  differ  by  tt/2  from  those  given  in  Alexander  &  Maddocks 
[1988].  The  scalar  pivot  or  steering  speed  is  defined  to  be  =  0,'(t).  The  vertical  component 
k  of  the  wheel  angular  velocity  is  0J(t)  =  +  0,  where  fi  is  the  angular  velocity  of  the  robot 

body  in  the  ambient  space. 

It  is  apparent  that  the  motion  of  a  simple  wheel  relative  to  the  body  of  the  WMR  is  math¬ 
ematically  characterized  by  the  two  scalars  pi  and  Si,  or  by  pi,  w,  and  an  initial  value  0j(O). 
Moreover,  the  steering  angle  and  rotation  rate  are  the  physically  relevant  controls.  The  two 
scalars  p,  and  0,  uniquely  characterize  the  axle  vector  a ,  in  the  robot  body.  The  whole  point  of 
this  paper  is  to  investigate  how  the  relative  motions  of  the  rolling  wheels  and  the  robot  body, 
determine  the  motion  of  the  robot  body  in  the  ambient  space. 

The  position  (or  configuration)  of  the  robot  body  is  determined  by  any  two  of  the  axle- 
point  space  coordinates  X^(f).  Alternatively  the  body  configuration  is  totally  specified  given  one 
coordinate  X,(t)  and  the  body  orientation  angle  ©&  between  the  reference  directions  fixed  in  the 
body  and  fixed  in  space.  The  angular  velocity  of  the  body  is  ft  =  0[,(t).  It  is  apparent  that 

©<(<)  =  ©6(f) +0<(t), 


(2.1) 


Alexander  <k  Maddocks:  Kinematics  and  control  of  WMRs  4 

and  by  differentiation, 

e'-(t)  =  n(t) +ci><(«).  (2.2) 

The  fact  that  the  robot  is  a  rigid  body  is  expressed  by  the  rigidity  conditions 

jt  |Xi(t)  -Xy(t)|  =  jt  |x,-  -xy|  =  0,  (2.3) 

If  the  axle-point  velocity  v,(t)  =  X'  (f)  is  introduced,  (2.3)  can  be  expressed  as 

(X,- (t)  -  Xy(f))  •  (v{(f)  -  vy (t))  =0,  i,j  =  1, . .  ,m.  (2.4) 


This  condition  can  be  reformulated  in  another  useful  way.  Let  otij  be  the  signed  angle  between 
the  reference  direction  in  the  body  and  the  line  joining  X;  to  Xy.  Consideration  of  components 
parallel  and  perpendicular  to  X<  —  Xy  reveals  that  (2.4)  can  be  rewritten: 

|vi  |  sin(0j  —  otij)  —  | v y  |  sin(0y  —  a<y)  =  0.  (2.5) 

One  of  the  most  elementary,  yet  elegant,  results  in  kinematics  (Chasles’  theorem  or  Descartes’ 
principle  of  instantaneous  motion)  asserts  that  conditions  (2.4)  imply  that  at  each  instant  the 
motion  of  a  planar  rigid  body  coincides  with  either  (i)  a  pure  rotation  about  some  point  (the 
instantaneous  center  of  rotation  or  ICR)  or  (ii)  a  pure  translation.  It  is  apparent  that  given 
n/0,  and  the  location  t  of  the  ICR  then  the  velocity  v<  of  any  point  X,-  is  determined  by 
v*  =  n(t  -  X.)  X  k.  Conversely,  any  two  unequal  velocities  v,  and  Vy  of  points  X,  and  Xy 
satisfying  (2.4)  determine  the  ICR  and  0.  When  0  =  0,  any  two  velocities  v;  and  Vy  are  equal, 
and  vice  versa. 

3.  Rolling 

The  condition  of  rolling  of  the  ith  wheel  relates  the  axle  vector  a,,  which  characterizes  the  wheel 
motion  relative  to  the  body,  to  the  axle-point  velocity  v<,  which  partially  describes  the  motion 
of  the  robot  body  in  space.  Explicitly,  the  rolling  conditions  are: 

(i)  the  directed  angle  from  the  axle-point  velocity  v,  to  the  axle  vector  a,-  is  ir/2, 

(ii)  |v,|  =  2 nRi  | a,- 1  =  2 irRiPi. 

In  mathematical  terms,  =  2 X  k.  Physically,  v<  lies  in  the  plane  of  the  wheel  and  its 
magnitude  is  determined  by  the  rotation  rate  of  the  wheel. 

In  some  mechanics  texts,  the  further  distinction  of  pure  rolling  is  made.  Pure  rolling  is 
said  to  occur  when  the  angular  velocity  Wj  of  the  wheel  has  no  vertical  component,  or  in  our 
notation,  +  Q  =  0.  This  last  concept  is  not  relevant  here,  and  the  term  pure  rolling  is  reserved 
to  emphasize  the  contrast  between  rolling  and  the  case  of  combined  rolling  and  slipping  that  is 
treated  in  §§7  and  8. 

The  condition  of  rolling  combines  naturally  with  the  notion  of  instantaneous  center  of  rota¬ 
tion.  If  the  instantaneous  motion  is  a  translation,  then  the  axle  vector  of  any  rolling  wheel  must 
be  orthogonal  to  the  translation  and  be  of  given  magnitude.  More  interestingly,  if  the  body  is 
rotating  with  instantaneous  velocity  Q  /  0,  and  the  ith  wheel  is  rolling,  then  the  ICR  must  be 
located  at  the  point 


X,-f27rf2jn  1a<. 


(3.1) 
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In  the  sequel  we  shall  simplify  notation  by  assuming,  with  no  loss  of  generality,  that  Ri  =  R}  = 
1/2tt.  This  assumption  merely  represents  a  scaling  of  each  axle  vector  a^.  Consequently,  any  two 
rolling  wheels,  say  the  ith  and  the  yth,  must  satisfy  the  rolling  compatibility  condition : 

fi(Xj  —  Xy)  =  fi(xi  —  xy)  —  (ay  —  at),  (3-2) 

which  expression  remains  valid  in  the  translational  case  0  =  0.  It  is  apparent  from  (3.2)  that 
once  axle  vectors  corresponding  to  two  rolling  wheels  are  given,  then  the  angular  velocity  0  and 
therefore  the  axle  vector  of  any  other  rolling  wheel  is  determined.  On  the  other  hand,  even  the 
first  two  axle  vectors  are  not  independent,  for  the  vector  equation  (3.2)  must  be  solvable  for  the 
scalar  0.  Note  that  the  condition  for  (3.2)  to  be  solvable  is  nothing  other  than  (2.4),  which, 
taken  with  the  rolling  conditions,  states  that  ay  —  a,-  must  be  parallel  to  x<  —  xy.  Provided  this 
solvability  condition  is  satisfied,  we  find  that 

n  _  (a»  ~  •  (*.'  ~  xy) 


(3.3) 


4.  Rigid  body  motion  and  controls 

The  development  of  the  previous  sections  demonstrate  that  if  an  initial  position  and  orientation 
of  a  WMR  are  known,  and  two  compatible  axle  vectors  are  given  as  functions  of  time,  then  the 
complete  motion  of  the  robot  is  determined.  This  fact  follows  from  (3.1)  and  (3.2),  which  give 
the  location  of  the  ICR  and  the  angular  velocity  in  terms  of  known  quantities  (x^  —  xy),  a,-,  ay 
and  the  initial  location.  Consequently,  the  location  can  be  updated,  and  the  cycle  repeated.  In 
this  section,  we  explicitly  find  the  differential  equations  implicit  in  this  argument.  The  variables 
that  will  be  adopted  are  the  space  coordinates  of  one  axle  point  X,  and  the  body  orientation 

0b.  There  are  four  associated  controls,  namely  two  steering  angles  and  two  rotation  rates.  These 

four  scalars  must  satisfy  the  solvability  condition  (2.5).  Although  only  three  of  the  four  scalars 
are  independent,  there  is  no  convenient  way  to  explicitly  eliminate  one  of  them. 

The  differential  equations  governing  the  motion  are 

X<=v<(  (4.1a) 

©(.  =  O-  (4.1b) 

However,  with  the  assumption  of  rolling,  the  space  coordinates  of  v<  can  be  expressed  in  terms 
of  6i  and  06;  also  (3.3)  provides  an  expression  for  fi.  Thus  (4.1)  can  be  recast  as  a  system  of 
three  first-order  equations  ((4.2a)  is  a  vector  equation  written  in  terms  of  two  coordinates): 

X-  =  p,(sin(0t-  +  ©(,),  cos(0j  -I-  ©6)),  (4.2a) 

®'b=  ~  [p<  cos(0i  -  a<y)  -  pj  cos(0y  -  a<y)] ,  (4.2b) 

where  x,y  =  |x,-  -Xy|.  Here  the  reference  direction  in  space  is  taken  to  lie  along  the  positive  Xx- 
axis.  Equations  (4.2)  taken  with  a  set  of  initial  conditions  determine  the  location  and  orientation 
of  the  robot  given  a  set  of  compatible  controls.  The  equation 


xy  =  Py(sin(0y  +  ©&),cos  (0y  +  0fc)), 


(4.3) 
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can  be  solved  for  Xy,  and  provided  that  (2.5)  is  satisfied,  |X;  (t)  —  Xy(t)|  is  automatically  con¬ 
stant.  Alternatively  Xy  can  be  constructed  from  the  formula 


Xy  —  X,'  +  Xij(coa{ctij  +  ©(>),  sin(a;y  +  ©b))-  (4.4) 

The  above  analysis  resolves  the  inverse  (or  control)  problem  en  passant.  For  if  X{  and  ©b 
are  prescribed  as  functions  of  time,  then  (4.2a)  can  be  regarded  as  two  coupled  transcendental 
equations  that  are  uniquely  solvable  for  the  controls  p,-  >  0,  and  &i  modulo  2n.  Moreover,  with 
given  X,-  and  ©b,  (4.4)  provides  an  expression  for  Xy.  Consequently,  (4.3)  provides  a  unique 
py  >  0,  and  6j,  modulo  2tt.  Similarly,  if  X,  and  Xy  satisfying  (2.2)  are  prescribed  as  functions  of 
time,  then  (4.4)  provides  ©b,  and  (4.2a)  and  (4.3)  can  be  solved  for  the  controls  as  before.  Any 
set  of  four  controls  obtained  in  this  way  are  automatically  compatible. 

More  specifically,  suppose  that  the  position  of  a  point  Y  of  the  WMR  and  ©b  are  prescribed 
as  functions  of  time.  There  is  no  loss  of  generality  in  assuming  that  Y  is  actually  an  axle  point 
Xt,  because  a  “virtual”  wheel  can  be  imagined.  From  (4.2a),  we  obtain 

Pi  =  |X'I,  (4.5a) 

0i= /3i-  ©b,  (4.56) 

where  Pi  is  the  known  angle  between  X'  and  the  Xa-axis  of  spatial  coordinates.  To  obtain  the 
controls  for  the  other  wheels,  differentiate  (4.4)  to  obtain 

Xy  =  X'  +  XijQ'b  (-  sin  (a<y  +  ©b),  cos(ayy  +  06)) , 
and  use  (4.2a)  and  (4.3)  to  obtain 

Pi  =  \Jp ?  +  *?y(©b)3  +  2p,x<y©;cos(^  +  ocij  +  2©,,). 

Thus  the  control  equations  are 

Pi  =  \/lX^l2  +  +  2|XJ|*jy©£  coB(ft  +  a,y  +  0b) 

/) .  _  *„„-i  {  lX(  |  sin  Pi  -  Xi}@'h  sin(a,y  +  ©6)  \  „ 

’  l|x;|cosA  +  xi)etco,K  +  0s)J'-8‘' 

Formulae  (4.6)  give  closed  form  expressions  for  the  controls  py,  0y  of  an  arbitrary  wheel  in  terms 
of  known  quantities.  In  particular,  there  are  no  differential  equations  involved. 

One  of  the  simplest  practical  problems  in  control  of  a  WMR  is  to  move  the  robot  body 
from  one  configuration,  as  represented  by  X,(0)  and  ©b(0)  say,  to  a  final  configuration  Xt(T) 
and  ©b(T).  The  above  analysis  demonstrates  that  to  each  rigid  body  trajectory  linking  the 
initial  and  final  configurations,  or  equivalently  to  each  pair  of  functions  X;(t)  and  06(t)  with 
correct  initial  and  final  values,  there  corresponds  a  unique  set  of  controls.  Consequently,  many 
algorithms  for  the  generation  of  controls  can  be  developed.  The  mathematically  simplest  choice 
involves  functions  X,-(f)  and  which  are  linear  in  time.  In  practice  other,  non-kinematic, 

considerations  enter  into  the  question  of  which  generating  algorithm  is,  in  some  sense,  best.  For 
example,  an  optimal  control  study  could  start  from  this  analysis,  and  would  further  determine 
which  set  of  controls  amongst  all  those  feasible  actually  minimized  some  functional,  perhaps  total 


(4.6a) 

(4.66) 
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work.*  Another  plausible  criterion  for  the  best  algorithm  might  be  the  most  robust  one,  i.e., 
the  control  for  which  errors  matter  the  least  (in  some  sense),  or  which  is  least  likely  to  lead  to 
slippage. 

The  analysis  of  this  section  takes  no  account  of  constraints  on  the  values  that  controls 
may  assume.  For  example,  most  practical  designs  will  involve  a  maximum  steering  lock,  and 
accordingly  the  control  9i  must  lie  in  some  predetermined  range.  The  development  presented  here 
certainly  allows  the  admissibility  of  a  given  trajectory  and  its  associated  controls  to  be  checked 
a  posteriori,  but  we  do  not  attempt  the  much  harder  problem  of  the  a  priori  characterization  of 
trajectories  that  are  attainable  with  a  given,  restricted  set  of  controls. 

5.  Practical  designs 

The  analysis  of  the  previous  sections  has  demonstrated  that  in  order  to  be  able  to  maneuver 
a  WMR  in  arbitrary  planar  rigid  body  motions,  it  is  both  necessary  (in  the  considered  class  of 
designs)  and  sufficient  that  there  be  two  driven  and  steered  simple  wheels.  At  first  sight  it  appears 
that  any  two-wheeled  robot  must  immediately  fall  because  it  makes  only  two-point  contact  with 
the  supporting  plane.  This  argument  is  invalid  because  the  center  of  gravity  of  the  robot  body 
may  lie  below  the  axle  attachment  points.  See,  for  example  Helmers  [1983a],  which  discusses 
a  robot  with  two  non-simple  wheels  that  are  mounted  in  non-vertical  planes.  Nevertheless,  the 
body  of  a  robot  with  only  two  wheels  will  be  susceptible  to  pendulum  like  oscillations  that  may 
be  excited  during  maneuvers.  Consequently,  we  consider  robots  with  wheel  configurations  that 
provide  an  intrinsically  more  stable,  three-point  contact  with  the  supporting  plane. 

When  the  third  support  point  is  provided  by  a  third  simple  wheel  and  rolling  occurs,  the 
analysis  of  §3  applies  to  completely  determine  the  third  axle  vector  in  terms  of  the  first  two 
axle  vectors.  Consequently,  of  the  six  controls  associated  with  three  simple  wheels  only  three 
are  independent.  If  there  is  an  error  in  the  activation  of  any  one  of  these  controls,  slippage 
must  occur.  Accordingly,  it  is  apparent  that  there  is  considerable,  essentially  redundant,  work 
required  to  calculate  compatible  controls.  There  are  three  plausible  approaches  to  the  removal 
of  this  redundancy. 

(i)  The  mechanics  of  the  maneuvering  robot  may  guarantee  that  some  of  the 
kinematic  compatibility  conditions  are  automatically  satisfied  when  some  of 
the  controls  are  left  passive.  For  example,  if  the  third  simple  wheel  is  steered 
correctly,  and  no  control  of  its  rotation  rate  is  made,  it  will  automatically  rotate 
at  the  correct  compatible  speed  in  response  to  the  frictional  forces  exerted  by 
the  surface.  In  contrast,  there  are  no  torques  that  would  automatically  orient 
the  third  wheel  and  allow  its  steering  control  to  be  left  passive.  Consequently, 
all  control  redundancies  cannot  be  eliminated  merely  by  leaving  redundant 
controls  passive. 

(u)  Some  mobility  of  the  robot  body  can  be  sacrificed  by  restriction  of  the  range 
of  some  of  the  controls.  Consider,  for  example,  a  two-wheeled  robot.  As  in  the 
argument  of  the  previous  paragraph,  the  forces  generated  by  one  non-vanishing 
rotation  rate  allows  the  other  rotation  control  to  be  left  passive.  Thus,  if  it  is 
accepted  that  pure  rotations  about  one  of  the  axle  points  (say  the  first)  cannot 
be  accessed,  then  only  three  controls  need  be  actuated.  However,  even  these 


*  Possibly  the  functional  minimized  is  the  total  work  of  the  designer. 
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three  controls  are  not  quite  independent;  whenever  0 2  =  c*i2,  equation  (2.5)  is 
an  active  constraint.  One  convenient  way  to  eliminate  even  this  compatibility 
condition  is  to  further  limit  the  accessible  class  of  rigid  body  motions  by  fixing 
the  steering  angle  of  the  second  wheel  with  9 2  ^  £*1.2;  i-  e.,  require  the  direction 
of  the  second  axle  vector  to  be  fixed,  and  not  to  be  parallel  to  the  line  between 
axle  points.  In  the  instance  of  two  wheels,  comparatively  few  motions  are 
excluded  by  fixing  one  axle.  However,  for  a  three-wheeled  robot,  the  steering 
control  of  two  wheels  would  have  to  be  held  constant  in  order  to  eliminate 
all  compatibility  conditions.  Having  two  fixed  wheels  limits  possible  rolling 
motions  to  the  unacceptably  small  set  of  circular  and  straight  line  trajectories, 
unless  the  fixed  axle  vectors  are  parallel  to  each  other,  and  to  the  line  between 
the  two  axle  points.  This  last  configuration  is  equivalent  to  having  two  simple 
wheels  mounted  on  the  same  axle  with  independent  rotation  rates  (such  as  the 
rear  wheels  of  an  automobile).  This  special  three  wheel  configuration  with  one 
steerable,  two  fixed,  one  driven,  and  two  passive,  wheels  is  a  viable  design  which 
was  adopted  for  the  WMR  Neptune  (Podnar,  Dowling  &i  Blackwell  [1984]).  Its 
kinematics  are  further  discussed  in  §6. 

(iii)  If  the  restricted  classes  of  planar  rigid  body  motions  described  in  (ii)  are  not 
acceptable,  then  a  mechanically  stable  robot  design  with  minimal  control  com¬ 
patibility  conditions  can  be  achieved  by  having  two  steered,  driven  wheels  and 
a  passive  castor.  For  our  purposes  a  castor  is  a  second  (small)  planar  rigid 
body  which  has  mounted  in  it  a  fixed  simple  wheel,  and  which  is  attached  to 
the  robot  body  at  a  free  pivot  point.  The  addition  of  a  castor  to  two  simple 
wheels  certainly  makes  the  robot  body  mechanically  stable.  The  four  controls 
associated  with  the  two  simple  wheels  must  satisfy  the  one  associated  com¬ 
patibility  condition,  but  the  castor  can  be  left  entirely  passive  with  no  rolling 
condition  being  violated.  This  last  fact  can  be  seen  from  a  modification  of  the 
arguments  of  the  previous  paragraph.  The  control  of  the  two  simple  wheels 
determines  the  motion  of  the  pivot  point  of  the  castor,  which  we  imagine  to 
be  the  axle  point  of  a  fictitious,  third,  simple  wheel.  Moreover,  the  axle  vector 
of  this  imaginary  wheel  is  determined.  But  the  imaginary  wheel  can  also  be 
regarded  as  being  mounted  in  the  castor  body,  and  as  such  it  completely  de¬ 
termines  the  motion  of  the  castor  body.  Furthermore,  provided  that  the  axle 
vector  of  the  wheel  fixed  in  the  castor  does  not  pass  through  the  imaginary 
wheel  (as  is  the  case  for  the  usual  design  of  castors),  the  arguments  of  (ii) 
above  apply  to  state  that  the  rotation  rate  of  the  castor  wheel  may  be  left 
as  a  passive  control.  Thus  the  castor  is  completely  passive.  Notice  that  the 
restrictions  on  possible  rigid  body  motions  that  arose  in  (ii)  imply  that  no 
matter  how  the  robot  body  is  maneuvered,  the  castor  will  not  undergo  certain 
rigid  body  motions.  However,  the  motion  of  the  WMR  itself  is  not  restricted, 
so  no  difficulty  arises.  It  is  apparent  that  passive  castors  are  irrelevant  to  the 

kinematics  of  a  WMR;  they  can  be  ignored  in  the  development  of  the  equations 
of  motion. 


All  of  the  above  described  configurations  involve  some  passive  controls.  It  should  be  re- 
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marked  that  additional  information  can  be  gained  if  these  passive  controls  are  actually  sensed. 
In  particular,  it  could,  in  principle,  be  verified  that  the  sensed  controls  are  compatible  with  pure 
rolling.  Such  sensing  provides  a  test  that  could  indicate  the  presence  of  slippage,  and  concomi¬ 
tantly  the  necessity  of  the  model  developed  in  §7.  The  analysis  required  to  exploit  sensed  data 
from  simple  wheels  is  straightforward,  but  the  problem  of  interpreting  sensed  data  from  a  castor 
is  more  complicated,  and  will  not  be  pursued  here. 

6.  Examples  of  rolling 

As  an  application  of  our  theory,  we  consider  the  elementary,  but  basic,  control  problem  of  moving 
a  WMR  from  one  prescribed  rigid  body  configuration  to  another.  We  first  consider  this  problem 
for  a  completely  mobile  robot,  that  is  for  a  robot  that  can  traverse  arbitrary  rigid  body  trajec¬ 
tories.  We  then  consider  the  problem  for  other  WMR  configurations  that  can  only  traverse  a 
restricted  class  of  rigid  body  trajectories. 

Consider  first  a  completely  mobile  robot.  Two  such  designs  are  the  CMU  Rover  (Moravec 
[1983])  which  has  three  driven,  steered  simple  wheels,  and  the  design  suggested  in  (5.iii),  namely 
a  WMR  with  two  driven,  steered  simple  wheels,  and  a  supporting,  passive  castor.  Because  the 
WMR  is  assumed  to  be  completely  mobile  there  is  no  question  of  which  trajectories  linking  the 
given  starting  and  ending  configurations  are  actually  accessible.  We  shall  therefore  pick  a  simple 
trajectory,  and  determine  the  controls  that  move  one  axle-point  X,  uniformly  along  the  X2- 
axii  from  (0,0)  to  (0,X)  in  unit  time,  whilst  simultaneously  rotating  the  robot  uniformly  from 
©6(0)  =  0°  to  ©fc(l)  =  0J.  We  set 


X,(t)  =  (0,0)  +  t(0,X),  (6.1a) 

©<>(*)  =  ©2(1  -t)  +  *©£.  (6.16) 

From  (4.5b),  64  =  &  —  06  =  -©°(1  —  t)  —  t&l-  Moreover,  from  (4.5a)  p,  =  |AT/|  =  X.  The 
controls  to  be  applied  at  another  simple  wheel  mounted  at  the  axle  point  Xy,  are  obtained  from 
equations  (4.6),  which  reduce  to 


*  =  \/ * 2  +  xUQl  -  ©b)2  +  2 Xxi}(el  -  eg)  cos(a,y  +  ©„(*», 


—  tan  1 


(  ®°)  sin(at-y  +  ©b(t))  ~\ 

VX+X,y(0l  -©O)cos(aiy  +  06(t))J  " 


(6.2a) 

(6.26) 


and  06(f)  is  given  by  (3.1b).  Here  x,y,  defined  in  §5,  and  a.y,  defined  in  §2,  are  constants 
determined  by  the  geometry  of  the  WMR.  If  any  of  the  simple  wheels  of  the  robot  are  not 
controlled  exactly  as  is  specified  by  (6.2)  slippage  of  the  wheels  must  occur,  and  the  analysis  of 
§§7  and  8  will  be  required. 

Another  practical  WMR  wheel  configuration  is  kinematically  equivalent  to  a  child’s  tricycle. 
The  kinematics  of  such  vehicles  is  considered  in  Alexander  &  Maddocks  [1988].  The  WMR 
Newt  (Hollis  [1977])  is  one  tricycle-like  WMR.  It  has  two  parallel  driven  wheels  that  are  fixed 
in  direction,  along  with  a  passive  castor  for  stability.  Another  existing  robot,  Neptune  (Podnar, 
Dowlmg  &  Blackwell  [1984]),  has  two  passive  wheels  fixed  in  direction,  and  a  parallel  driven, 
steered  wheel.  Because  tricycle-like  WMRs  have  a  fixed  axle,  they  are  not  completely  mobile’ 
and  therefore  arbitrary  rigid  body  trajectories  cannot  be  specified.  For  example,  such  a  WMR 
cannot  be  moved  parallel  to  the  fixed  axle.  However,  the  path  of  any  given  point  not  on  the  fixed 
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axle  can  be  prescribed  arbitrarily  (Alexander  &  Maddocks  [1988]).  In  particular,  a  WMR  with 
the  kinematics  of  a  tricycle  can  be  maneuvered  in  any  pure  rigid-body  rotation  where  the  center 
of  curvature  lies  on  the  extension  of  the  axle  vector. 

If  we  assume  there  is  no  further  restriction  on  a  tricycle- like  WMR,  such  as  maximum  steering 
lock,  it  is  possible  to  implement  a  simple  control  strategy  that  will  move  the  robot  between  two 
arbitrary  configurations.  The  strategy  we  implement  comprises  two  linked  rotations.  Thus 
the  WMR  travels  in  two  circular  arcs  with  constant  controls  on  each  arc.  Consider  any  circle 
through  the  fixed  axle  point  in  its  initial  position,  that  also  cuts  the  fixed  axle  vector  orthogonally. 
Consider  also,  the  corresponding  family  of  circles  through  the  final  position.  Of  the  family  of 
circles  through  the  final  position,  there  are  two  circles  which  are  tangent  to  the  initial  circle.  If  the 
WMR  is  steered  along  the  initial  circle  to  one  of  the  points  of  mutual  tangency,  and  then  steered 
along  the  second  circle,  the  robot  reaches  the  desired  final  position  in  the  correct  orientation. 
If  the  robot  is  steered  to  the  other  point  of  mutual  tangency,  and  then  along  the  corresponding 
second  circle,  the  final  position  reached  will  be  a  rotation  through  7r  of  the  desired  orientation. 

The  Stanford  Cart  (Moravec  [1980];  Moravec  [1983])  is  another  WMR  that  is  essentially 
kinematically  equivalent  to  a  tricycle.  It  has  two  parallel  driven  wheels  fixed  in  direction,  and  a 
pair  of  steerable  wheels,  configured  as  a  small  automobile.  For  our  purposes  the  only  difference 
between  the  Stanford  cart  and  a  tricycle  is  that  there  are  more  wheels,  and  therefore  more 
redundant  controls  that  must  be  accurately  activated  to  obviate  slippage.  For  WMRs  with  four 
or  more  wheels,  there  is  almost  certain  to  be  some  slippage  due  to  incompatibility  of  activated 
wheel  controls. 

7.  Slippage 

Sections  2-5  provide  an  analysis  that  uses  a  kinematic  model  of  pure  rolling  of  simple  wheels  to 
predict  the  motion  of  a  WMR.  The  analysis  does  not  explicitly  solve  Newton’s  laws  of  motion. 
Instead,  it  is  implicitly  assumed  that  the  undetermined  static  frictional  forces,  which  act  between 
the  wheels  and  the  supporting  plane,  are  sufficiently  large  as  to  be  able  to  balance  the  inertial 
forces,  which  arise  from  the  accelerations  of  the  WMR  predicted  by  the  rolling  model.  Naturally, 
the  validity  of  this  assumption  depends  upon  the  relative  sizes  of  the  inertial  and  static  frictional 
forces.  Whenever  the  inertial  forces  dominate,  the  rolling  model  ceases  to  be  valid.  The  subse¬ 
quent  motion  involves  a  phenomenon  we  call  skidding,  and  any  analysis  would  have  to  be  based 
on  the  full  system  of  Newton’s  laws.  On  the  other  hand,  the  rolling  model  also  fails,  even  if  the 
inertial  forces  are  negligible,  whenever  the  wheel  controls  are  incompatible,  i.  e.whenever  rolling 
conditions  (3.2)  are  not  satisfied.  The  subsequent  motion  involves  the  phenomenon  we  call  slip¬ 
page.  The  purpose  of  this  section  is  to  present  a  quasi-static  model  which  can  predict  motions 
that  include  slippage  due  to  wheel  incompatibility.  The  model  described  is  quasi-static  because, 
in  distinction  from  the  preceding  analysis,  the  inertial  forces  are  assumed  to  be  sufficiently  small 
as  to  be  totally  negligible. 

In  the  previously  described  model  of  rolling,  we  obtained  closed  form  expressions  for  the  axle 
point  velocities  v<  in  terms  of  the  controls,  as  represented  by  the  axle  vectors  a,-,  and  consequently 
obtained  differential  equations  that  describe  the  axle  point  trajectories.  Here  we  pursue  a  less 
ambitious  goal,  namely  to  find  an  algorithm  that  provides  the  axle  point  velocities  v,  given  the 
current  location,  as  represented  by  X;,  and  the  current  controls,  as  represented  by  a,'.  We  also 
consider  the  problem  in  which  one  or  more  drive  controls  are  left  passive.  In  this  case,  the 
direction  of  the  axle  vector  is  known,  but  its  magnitude  is  undetermined. 
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Our  first  step  is  the  construction  of  a  functional  that  represents  the  power  dissipated  by 
Coulomb  frictional  forces  during  any  planar  rigid  body  motion  of  the  WMR.  We  then  obtain 
a  quasi-static  evolutionary  system  for  the  WMR  by  positing  a  principle  of  least  power,  which 
states  that  the  actual  motion  of  the  WMR  is  the  one  minimizing  power  dissipation  at  all  times.  A 
closely  related  principle  was  introduced  by  Peshkin  and  Sanderson  (Peshkin  k  Sanderson  [1986]) 
to  describe  quasi-static  sliding  of  systems  of  particles,  but  rolling  was  not  considered  there. 
Our  approach  also  has  obvious  connections  with  Rayleigh’s  dissipation  function  (Lord  Rayleigh 
[1894-1896];  Whittaker  [1961,  p.  230]). 

Although  the  robot  body  is  no  longer  undergoing  pure  rolling,  it  is  still  moving  as  a  planar 
rigid  body,  and  therefore  has  an  angular  velocity  Cl.  It  should  be  recalled  that  either  0  =  0  and 
the  axle  point  velocities  V;  are  all  equal  (the  case  of  translation),  or  the  axle  point  velocities  are 
determined  by  v,  =  0(t  —  Xi)  X  k  where  X,  is  the  axle  point  location,  and  t  is  the  location  of 
the  instantaneous  center  of  rotation.  In  either  case,  we  can  write 

v^fT-nXiJxk,  (7.1) 


where  T  =  Ot  if  O  ^  0,  and,  when  O  =  0,  T  is  the  vector  satisfying  v,-  =  vy  =  T  x  k.  With  this 
definition,  T  behaves  smoothly  in  the  limit  0  — ♦  0. 

Consider  a  simple  wheel  that  is  undergoing  a  combination  of  rolling  and  slippage.  Define 
the  rolling  velocity  r;  by  r,  =  x  k  where,  as  before,  the  a,-  are  the  (scaled  with  a  factor  2a- Ri) 
axle  vectors,  and  k  is  the  vertical  unit  vector.  In  the  case  of  pure  rolling  r,  coincides  with  the 
velocity  v,  of  the  axle-point  X,  ,  but  in  the  case  of  combined  rolling  and  sliding  there  can  be  a 
discrepancy.  We  define  the  slippage  velocity  a(-  to  be  this  discrepancy.  Explicitly,  for  all  i, 

Si  =  v,  ~n  =  {(T-nXi)  -ai}  xk,  (7.2) 

It  should  also  be  remarked  that  the  wheel  rotation  rate  pi  satisfies 

*  =  M=|a<|.  (7.3) 

According  to  Coulomb’s  law,  the  sliding  friction  fj  between  a  particle  and  a  surface  is  given 
by  the  expression 


where  Si  is  the  relative  velocity  of  the  particle  and  the  surface,  and  a{  >  0  is  the  product  of 
the  normal  load  and  the  coefficient  of  friction.  Then  the  total  instantaneous  power  P  being 
dissipated  by  sliding  friction  is  the  sum  of  the  scalar  products 

-P=-5If*'8i>  (7.5) 

* 

or,  by  (7.4), 


*-X>W-  (7-6) 

t 

Of  course,  the  construction  of  (7.6)  applies  to  the  contact  points  between  the  wheels  of 
the  WMR  and  the  supporting  plane,  with  the  relative  velocity  being  the  slippage  velocity  s,-. 
Consequently,  (7.2)  can  be  used  to  rewrite  (7.6)  as 

P  =  P(n,T;Xi,a,)  =  ^ai|(T  -  DXi  -  a,-)  x  k 


(7.7) 
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We  call  (7.7)  the  dissipation  or  friction  functional.  Because  k  is  a  unit  vector,  and  all  other 
vectors  appearing  are  perpendicular  to  k,  the  dissipation  functional  can  also  be  written  in  the 
form: 

P  =  P(  n,T;Xi,a,)  =  ^«<|(T— nXi-a^l.  (7.8) 

t 

The  quasi-static  dynamical  system  governing  the  WMR  motion  with  slippage  is  obtained 
by  taking  the  values  of  fi  and  T  that  realize  the  minimum  of  the  slippage  functional.  When 
fl^O,  knowledge  of  fi  and  T  provides  the  angular  velocity  and  instantaneous  center  of  rotation 
t.  When  fi  =  0,  knowledge  of  T  provides  the  common  translational  velocity  v.  During  the 
minimization,  the  quantities  X;  and  a;  are  assumed  known  and  are  held  fixed.  In  the  case  of 
a  passive  drive  control  the  interpretation  is  slightly  different.  Then  the  axle  vector  a j  is  of  the 
form  p,Aj,  where  A,-  is  a  given  unit  vector,  but  the  rotation  rate  Pi  is  a  priori  unknown.  Any 
unknown  rotation  rates  are  regarded  as  additional  variables  over  which  the  slippage  functional 
must  also  be  minimized. 

The  dissipation  functional  is  not  smooth  whenever  T  =  fiXj  +a<.  These  cases  arise  precisely 
when  the  ith  wheel  is  undergoing  pure  rolling.  The  lack  of  smoothness  is  explained  in  physical 
terms  below.  When  the  dissipation  functional  is  smooth,  the  first  order  necessary  conditions  for 
local  minimization  comprise  equating  the  partial  derivatives  of  (7.7),  or  equivalently  (7.8),  with 
respect  to  T  to  fi  equal  to  zero.  The  resulting  slippage  evolution  equations  are 


E(T  -  nx<  -  a<)  x  k  „  ^  T  -  fix,  -  &i 
,  |  (T  —  fiX<  —  a<)  X  k|  ~  °  ~  2-j  “‘"iT-fiXi  -  a,-|  ’ 


and 


T  —  fiX,-  —  a,- 
|T  —  fiXi  —  a,- 1 


(7.9) 


(7.10) 


Equations  (7.9)  and  (7.10)  possess  natural  physical  interpretations.  The  vector  equation 
(7.9)  expresses  the  balance  of  frictional  forces  on  the  robot:  from  (7.2)  and  (7.4), 


Ef*  =  °- 


(7.11) 


The  scalar  equation  (7.10)  expresses  the  balance  of  the  moments  of  all  frictional  forces  acting  on 
the  body: 

J3(X  -  X,)  x  fi  =  0  (7.12). 

t 

where  X  is  an  arbitrary  point.  Use  of  (7.2),  (7.9)  and  the  expansion  formula  for  vector  triple 
products  reduces  the  moment  balance  (7.12)  to  condition  (7.10). 

The  physical  interpretation  of  (7.9)  and  (7.10)  rigorously  justifies,  in  the  smooth  case,  the 
equivalence  of  the  principle  of  minimal  dissipation  and  the  usual  approach  to  quasi-static  evo¬ 
lutionary  systems,  namely  a  balance  of  all  forces  and  moments  except  inertial  terms.  We  also 
take  the  equivalence  of  smooth  cases  to  be  evidence  supporting  the  validity  of  the  principle  of 
least  dissipation  when  the  friction  functional  is  not  smooth.  There  are  actually  sound  physical 
reasons  to  expect  nonsmoothness  of  the  friction  functional.  As  was  pointed  out  previously,  nons¬ 
moothness  arises  precisely  when  one  or  more  wheels  are  undergoing  pure  rolling.  For  pure  rolling 
the  possible  frictional  forces  between  the  wheel  and  the  supporting  plane  are  restricted  only  by 
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a  static  friction  law.  In  particular,  the  direction  and  magnitude  of  the  frictional  forces  are  not 
completely  determined,  and  an  explicit  force  and  moment  balance  cannot  be  expected. 

When  passive  controls  occur,  so  that  a  wheel  is  steered  but  not  driven,  further  first  order 
conditions  are  obtained  by  our  equating  to  zero  the  partial  derivatives  of  the  slippage  functional 
with  respect  to  the  passive  rotation  rates.  More  precisely,  for  each  passive  wheel  i, 

a<((T  -  fiXi  -  p<Ai)  x  k)  •  (A i  x  k)  =  0,  (7.13) 

However  this  equation  is  equivalent  to  the  geometrical  condition  that  any  passive  rotation  rate 
must  be  chosen  such  that  the  slippage  velocity  is  orthogonal  to  the  rolling  velocity,  or  s <  •  U  =  0. 
Hence  there  is  a  minimization  principle  at  each  passive  wheel;  for  given  V;  and  a,-,  the  slippage 
vector  Si  has  minimal  length. 


8.  Examples  involving  slippage 

Before  a  general  analysis  of  the  kinematics  of  WMRs  with  slippage  can  be  attempted,  the  dis¬ 
sipation  functional  must  be  better  understood.  We  here  present,  as  an  illustrative  example,  an 
analysis  of  the  case  of  two  wheels. 

As  was  noted  in  §5,  when  a  WMR  has  two  wheels  both  of  which  are  driven  and  steered, 
then  there  is  one  constraint  (3.2)  on  the  controls  that  must  be  satisfied  in  order  that  rolling 
is  possible;  if  this  constraint  is  violated,  slippage  must  occur.  We  shall  analyze  the  motion, 
including  slippage,  for  such  two  wheeled  robots.  The  remarks  concerning  castors  made  in  §5 
apply  with  equal  force  here  to  imply  that  castors  have  no  effect  on  the  kinematics  of  a  WMR. 
Thus  it  is  apparent  that  the  subsequent  analysis  also  applies  to  robots  with  two  simple  wheels 
and  passive  stabilizing  castors. 

Denote  the  wheels  by  i  =  1,  2.  Recall  that  the  a<  >  0  are  the  product  of  the  normal  load  on 
the  ith  wheel  and  the  coefficient  of  friction.  There  are  two  cases  to  consider:  (i)  the  symmetric 
case  ai  =  <*2  and  (ii)  the  asymmetric  case  a\  ^  a2.  In  the  case  of  two  wheels,  (7.9)  becomes 
the  weighted  sum  of  two  unit  vectors,  and  accordingly  has  a  solution  if  and  only  if  the  two  unit 
vectors  are  anti-parallel,  and  Qi  =  ot  2-  Hence,  in  the  asymmetric  case,  the  slippage  evolution 
equations  have  no  solution.  Therefore  any  minimum  must  occur  where  the  dissipation  functional 
is  not  differentiable  and  (7.8)  must  be  investigated  directly. 

Consider  the  asymmetric  case  first.  Suppose  that  c*i  <  a 2.  Here  the  functional  (7.8)  is 
minimized  when 

T  =  f2X2  +  a,',  (®-l) 

and  (7.8)  reduces  to 

ai|n(Xi  —  X2)  +  (ai  —  a2)|.  (8-2) 

From  this  last  expression  we  determine  that 


(ai-a2)  (Xx  -X2) 

|Xi  —  x2|2 


(8.3) 


Remarkably,  the  angular  velocity  17  in  (8.3)  that  is  predicted  by  the  slippage  model  coincides 
with  the  angular  velocity  of  (3.3)  predicted  by  the  rolling  model.  Notice  that  the  second  wheel 
undergoes  pure  rolling  in  such  a  manner  that  the  slippage  of  the  first  wheel  is  minimized.  When 
(ai  —  a2)  •  (Xx  —  X2)  =  0,  then  12  =  0,  and  the  motion  is  a  translation.  Furthermore,  the 
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minimal  dissipation  is  strictly  positive,  unless  ai  —  &2  is  parallel  to  Xi  —  X2,  in  which  case  the 
compatibility  condition  for  pure  rolling  for  both  is  satisfied,  P  =  0,  and  no  slippage  occurs. 
Consider  now  the  symmetric  case.  Here  (7.9)  is  satisfied  whenever 

T  =  (l  —  7)(flXi  +  ai)  +  7(0X2  +  82),  (8-4) 


for  0  <  7  <  1;  that  is,  whenever  T  lies  on  the  straight  line  segment  between  OXi  +  ai  and 
OX2  +  a2-  Using  (7.10),  we  determine  that 


(a1-a2)(X1-X2) 

|X!-X2|2 


(8.5) 


as  before.  Note  that  although  O  is  completely  determined,  condition  (8.4)  does  not  provide  a 
unique  solution  for  T.  Consequently,  neither  the  instantaneous  center  of  rotation  t  if  0  ^  0, 
nor  the  translational  velocity  v  if  0  =  0,  is  uniquely  determined.  In  practice,  it  might  be 
expected  that  the  load  on  one  wheel  or  the  other  will  dominate  at  any  particular  time  so 
that  that  wheel  will  experience  pure  rolling  at  that  instant.  Small  perturbations  could  cause  this 
dominance  to  jump  intermittently  from  one  wheel  to  the  other,  and  the  ICR  of  the  WMR  will 
vary  discontinuously.  * 
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